Decoherence by a spin thermal bath: Role of the spin-spin interactions and initial 

state of the bath 
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We study the decoherence of two coupled spins that interact with a spin-bath environment. It 
is shown that the connectivity and the coupling strength between the spins in the environment 
are of crucial importance for the decoherence of the central system. For the anisotropic spin-bath, 
changing the connectivity or coupling strenghts changes the decoherence of the central system from 
Gaussian to exponential decay law. The initial state of the environment is shown to affect the 
decoherence process in a qualitatively significant manner. 

PACS numbers: 03.67.Mn 05.45.Pq 75.10.Nr 



I. INTRODUCTION 

Understanding the decoherence in quantum spin sys- 
tems is a subject of numerous works (for reviews, see 
Refsii^). The issue seems to be very complicated and 
despite many efforts, even some basic questions about 
character of the decoherence process are unsolved yet. 
Due to the interactions with and between the spin of 
the bath, an analytical treatment can be carried out in 
exceptional cases, even if the central systems contains 
one spin only. Recent work suggests that the internal 
dynamics of the environment can be crucial to the de- 
coherence of the central system 3 i 4 i 5 ' 6 i 7 i 8 i 9 i 10 ' n i 12 i 13 ' 14 i 15 . 
In this paper, we present results of extensive simulation 
work of a two-spin system interacting with a spin-bath 
environment and show that the decoherence of the two- 
spin system can exhibit different behavior, depending on 
the characteristics of the coupling with the environment, 
the internal dynamics and the initial state of the latter. 
We also provide a simple physical picture to understand 
this behavior. 

In general, the behavior of an open quantum system 
crucially depends on the ratio of typical energy differ- 
ences of the central system SE C and the energy E ce which 
characterizes the interaction of the central system with 
the environment. The case SE C <C E ce has been stud- 
ied extensively in relation to the "Schrodinger cat" prob- 
lem and the physics is quite clear 16 i 17 : As a result of 
time evolution, the central system passes to one of the 
"pointer states"— which, in this case, are the eigenstates 
of the interaction Hamiltonian H ce . The opposite case, 
5E C ^> E ce is less well understood. There is a conjecture 
that in this case the pointer states should be eigenstates 
of the Hamiltonian H c of the central system but this 
has been proven for a very simple model onljsi^. On the 
other hand, this case is of primary interest if, say, the 
central system consists of electron spins whereas the en- 
vironment are nuclear spins, for instance if one considers 
the possibility of quantum computation using molecular 
magnet a 19 ' 20 . 



II. MODEL 

We consider a generic quantum spin model described 
by the Hamiltonian H = H c + H ce + H e where H c — 
— JSi • S2 is the Hamiltonian of the central system and 
the Hamiltonians of the environment and the interaction 
of the central system with the environment are given by 

N-l N 

H * = -E E E«#W 

i—l 1 a 

2 N 

i—l j—1 a 

respectively. The exchange integrals J and f^™- de- 
termine the strength of the interaction between spins 
S„ = (S®, S%, S*) of the central system, and the spins 
I" = (^n^ni^n) m the environment, respectively. The 
exchange integrals control the interaction of the cen- 
tral system with its environment. In Eq. ([1]), the sum 
over a runs over the x, y and z components of spin- 
1/2 operators S and I. In the sequel, we will use the 
term "Heisenberg-like" H e (H ce ) to indicate that each 
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(A^) is a uniform random number in the range 



% <3 

[— f2|,f2] ([—A, A]), Q, and A being free parameters. In 
earlier wor k 14 ' 15 , we found that a Heisenberg-like H e can 
induce close to perfect decoherence of the central system 
and therefore, we will focus on this case only. 

The bath is further characterized by the number of en- 
vironment spins K with which a spin in the environment 
interacts. If K = 0, each spin in the environment in- 
teracts with the central system only. K = 2, K = 4 or 
K = 6 correspond to environments in which the spins are 
placed on a ring, square or triangular lattice, respectively 
and interact with nearest-neighbors only. If K — N — 1 , 
each spin in the environment interacts with all the other 
spins in the environment and, to give this case a name, 
we will refer to this case as "spin glass" . 

If the Hamiltonian of the central system H c is a pertur- 
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bation relative to the interaction Hamiltonian H ce , the 
pointer states are eigenstates of H C J^. In the opposite 
case, that is the regime |A| <C \ J\ that we explore in this 
paper, the pointer states are conjectured to be eigenstates 
of H^.. The latter are given by |1) = \Ti) = |TT), |2) = 
\S) = (|U) - Iit»/V^, |3) = \T Q ) = (ITI) + |IT))/V2, 
and |4) = |T_ X ) = |||), satisfying H C \S) = (3J/4)|5) 
and H c \T t ) = (-J/4)\T t ) for i = -1,0, 1. 

The simulation procedure is as follows. We generate 
a random superposition |</>) of all the basis states of the 
environment. This state corresponds to the equilibrium 
density matrix of the environment at infinite tempera- 
ture. Alternatively, to study the effect of the thermal 
state of the environment on the decoherence processes, 
we take the 2 state of the environment to be its ground 
state. The spin-up - spin-down state (|T0) is taken as 
the initial state of the central system. Thus, the initial 
state of the whole system reads \\fr(t = 0)) = |||) \<j>) and 
is a product state of the state of the central system and 
the initial state of the environment which, in general is 
a (very complicated) linear combination of the 2 N basis 
states of the environment. In our simulations we take 
N = 16 which, from earlier wor k 14 i 15 , is sufficiently large 
for the environment to behave as a "large" system. 

For a given, fixed set of model parameters, the time 
evolution of the whole system is obtained by solving the 
time-dependent Schrodinger equation for the many-body 
wave function |\P(t)), describing the central system plus 
the environment 21 . It conserves the energy of the whole 
system to machine precision. We monitor the effects of 
the decoherence by computing the the matrix elements 
of the reduced density matrix p (t) of the central system. 

As explained earlier, in the regime of interest |A| <C 
\J\, the pointer states are expected to be the eigenstates 
of the central systems. Hence we compute the matrix el- 
ements of the density matrix in the basis of eigenvectors 
of the central system. We also compute the time depen- 
dence of quadratic entropy S c {t) = 1 — Trp 2 (t) and the 
Loschmidt echo L (t) = Tr (p (t) po (i))^, where po (t) is 
the density matrix for H ce = 0. 

III. ISOTROPIC COUPLING TO THE BATH 

If the interaction between the central system and en- 
vironment is isotropic we have [H c , H ce ] = 0. Then, as 
shown in the Appendix, the expressions of the reduced 
density matrix p (t) and the Loschmidt echo L (t) sim- 
plify. Indeed, if Ag = A 2 ( J = Ag = A for all 
then 

N 

ff ce = -A(S 1 + S 2 )-^I j (2) 

i=i 

commutes with H c and it follows that the decoherence 
process of the central system is determined by H ce , H e , 
the initial state of whole system |\&(to)), and the eigen- 
states of the central system (see Eq. (| L5[) and (|TB| in 
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FIG. 1: (Color online) The time evolution of the real part of 
the off-diagonal element pn (left panel) and the diagonal ele- 
ments pn, . . . , /944 (right panel) of the reduced density matrix 
of a central system (with J = —5), coupled via an isotropic 
Heisenberg interaction H ce (A = —0.075 ) to a Heisenberg- 
like environment H e (Q — 0.15) with different connectivity: 
(a) K = 0; (b) K = 2; (c) K = 4; (d) K = 6; (e) K = N - 1. 

the Appendix). In other words, in this case, L (t) and 
\p{t)\ do not dependent on the J, the interaction be- 
tween the spins in the central system. Furthermore, if 
we take the interactions between the environment spins 
to be isotropic, that is, — — fl\ z j = flij for all 
then 

N-l N 

II, X E ■ Ii (3) 

commutes with H cei and therefore H e has no effect on 
the decoherence process (see Eq. ([T5]l in the Appendix). 

In Fig. Q] and Fig. [21 we show the time evolution of the 
elements of the reduced density matrix p (t) for different 
connectivity K and f2, for the case that H ce is an isotropic 
Heisenberg model, i.e., A^ = A^ = Aj*? = A for all 
i,3- 
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FIG. 2: (Color online) The time evolution of the off-diagonal element P23 of the reduced density matrix of a central system (with 
J = —5), interacting with a Heisenberg-like environment H e via an isotropic Heisenberg Hamiltonian H ce (with A = —0.075 ) 
for the same geometric structures in the environment: (a,b) K = 2 and (c,d) K = TV — 1. The number next to each curve is 
the corresponding value of fi. 



If |A| ^> \/~KVL, in agreement with earlier wor k 23 ' 24 , 
we find that in the absence of interactions between the 
environment spins (yKQ = 0) and after the initial de- 
cay, the central system exhibits long-time oscillations (see 
Fig. da) (left)). In this case and in the limit of a large 
environment, we have 2 ^ 



Re p 2 3 (*) = 



l-bt 2 



cos cut, 



(4) 



where b = 7VA 2 /4, c = 6/2 and u = J - A. Equa- 
tion (HJ) clearly shows the two-step process, that is, after 
the initial Gaussian decay of the amplitude of the oscil- 
lations, the oscillations revive and their amplitude levels 
o&^. Due to conservation laws, this behavior does not 
change if the environment consists of an isotropic Heisen- 
berg system = f2 for all a, i and j), independent of 
K. If, as in Ref. 23 , we take A\ x ) = A^ = A$ € [0, A] 



random instead of the identical, the amplitude of the 
long-living oscillations is no longer constant but decays 
very slowly^ 3 - (results not shown). 

If |A| w VaA, the presence of Heisenberg-like inter- 
actions between the spins of the environment has little 
effect on the initial Gaussian decay of the central sys- 
tem, but it leads to a reduction and to a decay of the 
amplitude of the long-living oscillations. The larger K 
(see Fig. QIb-e)(left)) or ft (see Fig. H£a,c)), the faster 
the decay is. Note that for the sake of clarity, we have 
suppressed the fast oscillations by plotting instead of the 



real part, the absolute value of the matrix elements. 

If |A| <C \fKQ, keeping K fixed and increasing 
smoothly changes the initial decay from Gaussian (fast) 
to exponential (slow), and the long-living oscillations are 
completely suppressed (see Fig. HJb,d)). For large fi, the 
simulation data fits very well to 



\P23 (t)\ 



1 



-A K (n)t 



(5) 



with A K (fi) w VLA Kl A 2 = 9.13 and A N - X = 26.73. 
Note that, in principle, a closed quantum system cannot 
exhibit exponential decay^. The fact that we observe a 
decay that is well described by a single exponential may 
be the result of tracing out the degrees of freedom of 
an environment which initially is in a state of random 
superposition of the basis states. 

Physically, the observed behavior can be understood as 
follows. If | A| « ^/KQ, a bath spin is affected by roughly 
the same amount by the motion of both the other bath 
spins and by the two central spins. Therefore, each bath 
spin has enough freedom to follow the original dynamics, 
much as if there were no coupling between bath spins. 
This explains why the initial Gaussian decay is insensitive 
to the values of K or il. After the initial decay, the 
whole system is expected to reach an stationary state, but 
because of the presence of Heisenberg-like interactions 
between the bath spins, a new stationary state of the bath 
is established, suppressing the long-living oscillations. 

For increasing K, the distance between two bath spins, 
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FIG. 3: (Color online) Same as Fig. [2]except that H ce is Heisenberg-like and A = 0.15. 



denned as the minimum number of bonds connecting the 
two spins, becomes smaller. For instance, for K = 2, 
this distance is (N — 2) /2, and for K = N — 1, it is zero. 
Therefore, for fixed £1 and increasing K the fluctuations 
in the spin bath can propagate faster and the evolution to 
the stationary state will be faster. Similarly, for fixed K , 
increasing the coupling strength between the bath spins 
will speed up the dynamics of the bath, that is, the larger 
O the faster will be the evolution to the stationary state. 

In the opposite case |A| <C y/KQ, H ce is a small per- 
turbation relative to H e and the coupling between bath 
spins is the dominant factor in determining the dynamics 
of the bath spins. Therefore, by increasing K or fi, the 
bath spin will have less freedom to follow the dynamics 
induced by the coupling to the two central spins, the in- 
fluence of the bath on the central system will decrease, 
and the (exponential) decay will become slower. 

According to the general picture of decoherence^, for 
an environment with nontrivial internal dynamics that 
initially is in a random superposition of all its eigen- 
states, we expect that the central system will evolve to 
a stable mixture of its eigenstates. In other words, the 
decoherence will cause all the off-diagonal elements of 
the reduced density matrix to vanish with time. In the 
case of an isotropic Heisenberg coupling between the cen- 
tral system and the environment, H c commutes with the 
Hamiltonian H , hence the energy of the central system is 
a conserved quantity. Therefore, the weight of the singlet 
\S) in the mixed state should be a constant (1/2), and 
the weights of the degenerate eigenstates |T ), |T_i) and 
\Ti) are expected to become the same (1/6). As shown 
in Fig. [TJb-e) (right), our simulations confirm that this 



picture is correct in all respects. 

IV. ANISOTROPIC COUPLING TO THE BATH 

In order to clarify the role of K and fl, we change the 
coupling between the central system and the bath from 
Heisenberg to Heisenberg-like. From a comparison of the 
data in Fig. [5] and Fig. [31 it is clear that the roles of K 
and Q are the same in both cases, no matter whether the 
coupling to the bath is isotropic or anisotropic. How- 
ever, there are some differences in the decoherence pro- 
cess. The most important parameter determining the 
decoherence process is the ratio of the typical interaction 
energy A to the mean-square energy of interactions in 
the the thermal bath, y/~KQ. 

If |A| y/KQ : in the presence of anisotropic inter- 
actions between the central system and the environment 
spins, even in the absence of interactions between the 
bath spins, the second step of the oscillations decays and 
finally disappear as K increases. This is because the 
anisotropic interactions break the rotational symmetry 
of the coupling between central system and environment 
which is required for the long-living oscillations to per- 
sist. 

If |A| < y/Ktt, |p23 (t) I can still be described by 
Eq. ([5]) , but now Ak (fi) is no longer a linear function of 
O. For anisotropic H cei the energy of the central system 
is no longer a conserved quantity. Therefore there will 
be energy transfer between the central system and the 
environment and the weight of each pointer state (eigen- 
state) in the final stable mixture need not be the same 
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FIG. 4: (Color online) The time evolution of the the entropy 
S c (t) and Loschmidt echo L (t) of a central system (with 
J = —5), interacting with a Heisenberg-like environment H e 
(with different fi) via a Heisenberg (a,b, A = —0.075) or 
Heisenberg-like (c,d, A = 0.15) Hamiltonian H C e for the case 
K = 2. The number next to each curve is the corresponding 
value of Q. 

for all K or f2. 

For a change, we illustrate this point by considering 
the quadratic entropy S c (t) and Loschmidt echo L(t). 
We expect that these quantities will also dependent of 
the symmetry of the coupling between central system 
and the spin bath. In Fig. [51 we present results for 
large and K = 2, confirming this expectation. For 
isotropic (Heisenberg) H ce and perfect decoherence (zero 
off-diagonal terms in the reduced density matrix) we ex- 
pect that max. t S c (t) = 1- [(l/2) 2 + 3x (1/6) 2 ] = 2/3, in 
concert with the data of Fig. Ufa)). For Heisenberg-like 
H ce , max( S c (t) will depend on the coupling strengths 
and as shown in Fig. [He), we find that maxtS c (t) — 
1 — 4 x (1/4) 2 = 3/4, corresponding to the case that all 
the diagonal elements in the reduced density matrix are 
the same (1/4) and all other elements are zero. 

V. DISCUSSION AND CONCLUSIONS 

In the foregoing, we have compared y/KQ to |A| to 
distinguish different regimes. As a matter of fact, yKQ 
does not completely characterize the decoherence pro- 
cess, but it can be used to characterize its time scale. In- 
deed, as shown in Fig. for different \f~K and fl but the 
same value of \AKT2, the the time evolution of L(t) is very 
similar. Note that if y/~KQ, increases (compare Fig. [SJi to 
Fig. [5ji), the differences between the Loschmidt echoes 
increase. Additional simulations (results not shown) in- 
dicate that this differences are fluctuations that are due 
to the particular realization random parameters used in 
the simulation. 

In conclusion, for a spin-bath environment that ini- 



FIG. 5: (Color online) The time evolution of the Loschmidt 
echo L(t) of a central system (with J = —5), interact- 
ing with a Heisenberg-like environment H e via a Heisenberg 
(A = —0.075) Hamiltonian H ce . In each panel, the val- 
ues of y/~KQ are the same: (a) vKQ, = 0.1-y/N — 1, (b) 
•/KQ. = 0.15VA- 1, (c) VWn = 0.25VA- 1, and (d) 
y/KQ ee y/N — 1. The different lines in each pannel corre- 
spond to different K. Solid (black) line: K — 2; Dashed (red) 
line: K — 4; Dotted (green) line: K = 6, and dash-dotted 
(blue) line: K = N - 1. 



tially is in a random superposition of its basis states, we 
have shown how a pure quantum state of the central spin 
system evolves into a mixed state, and that if the interac- 
tion between the central system and environment is much 
smaller than the coupling between the spins in the cen- 
tral system, the pointer states are the eigenstates of the 
central system. Both these observations are in concert 
with the general picture of decoherenceii. Furthermore, 
we have demonstrated that, in the case that the environ- 
ment is a spin system, the details of this spin system are 
important for the decoherence of the central system. In 
particular, we have shown that for the anisotropic spin- 
bath, changing the internal dynamics of the environment 
(geometric structure or exchange couplings) may change 
the decoherence of the central spin system from Gaussian 
to exponential decay. 

Finally, we would like to compare the present results 
with those of our earlier work in which we focussed on the 
case in which the environment is initially in its ground 
state and demonstrated that, apart from the strength 
of different interactions, also their symmetry and the 
amount of entanglement of the ground state of the cen- 
tral system affects the decoherenc o 14 i 15 . To facilitate the 
comparison, in Fig. [5] we present some new data of the 
Loschmidt echoes for different K but for fixed y/KCl. 
Comparison of Fig. [5] with Fig. [5] indicates that if the 
environment is initially in its ground state, the decoher- 
ence process is qualitatively different from the one ob- 
served in the case that the initial state of the environ- 
ment is a random superposition. Roughly speaking, it 
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of the whole system e can be represented as 



FIG. 6: (Color online) The time evolution of the Loschmidt 
echo Lit) of a central system (with J = —5), interacting 
with a Heisenberg-like environment H e via a Heisenberg (A = 
—0.075) Hamiltonian H ce . The environment spins are initially 
prepared in the ground state. The different curves correspond 
to different K, but ^KQ. = 0.15VN - 1 is fixed. Solid (black) 
line: K — 2; Dashed (red) line: K = 4; Dotted (green) line: 
K — 6, and dash-dotted (blue) line: K = N - 1. 



is more difficult for the central system to change from 
a pure quantum state to a classical, mixed state, which 
is of course consistent with the fact that the quantum 
effects become more prominent as the temperature de- 
creases. In particular, from Fig. [5] it is clear that "JKQ, 
is not enough to characterize the qualitative behavior of 
the Loschmidt echo for the cases shown. 

The difference between the cases of an environment at 
low-temperatur o 14 ' 15 and a high-temperature (chaotic) 
environment considered in the present paper is most im- 
portant for the systems with very large connectivity. In 
the latter case, the ground state of the environment is 
a quantum spin-glass which is a very effective source 
of decoherenc o 14 ' 15 . At the same time, for the case of 
infinite temperature of the bath considered in this pa- 
per, this case is not very special when compared to the 
case of short-range interactions within the environment 
(see Fig. [SJ. It would be of interest to see if, as the 
temperature decreases, the decoherence process changes 
as the environment goes into the spin-glass state (at 
T oc y/KQ,), a problem that we leave for future research. 



VI. APPENDIX 



Consider a generic quantum model described by the 
Hamiltonian H = H c + H ce + H e , where H c and H e 
describe the central system and the bath respectively 
([H c ,H e ] = 0), and H ce describes the coupling between 
them. If [H c , H ce ] = 0, then the time evolution operator 



-iHt -iH c t -i(H C e+H c )t 



(6) 



Denote the eigenstates and corresponding eigenvalues of 
the central system by {\k}} and {Ek}, that is, H c \k) — 
Ek \k). The initial state (\f{to))) of the central system 
can be represented as \ip(t )) = J2k a k |&)- For an iso- 
lated central system (H ce = 0), the time evolution of the 
density matrix of the central system is given by 



PO (t)=Y, e ~ l{Ek ~ El)ta >< a l\ k ) 
k,l 



(7) 



If the central system is coupled to the bath Q(j) (to))), the 
initial state of the whole system can be represent as 



|*(io)) =$>*!*> I* (<o)>, 



(8) 



and the state at later time t is 
Mi))) = e- iHt \V(to)) 

= J2 e ~ iEkta ke- i{HcM)t \k)\d>(t )). (9) 
k 

As [H c ,H ce ] = 0, we have H ce \k) \(j> (t )) = 
\k) M k |</> (to)), therefore 

e -i(»~+ff.)*|A;)^(fo)) 

= E ( " rtr(g T + Her i fe )i^( f °)) 

ml 

m 

* — 4 ml 

m 

= \k)e-< M *+ H ^\<f>{t )) 

= |fc> 1^* (*)> , (10) 
where we introduced 

\M))=e-* M * +H °» |0 (t )). (if) 
Hence, the state at time t becomes 

|*(i)) = 5> fc e-^ 4 |fc) |0 fe (t)) . (12) 
k 

The density matrix p (t) of the whole system is 

P (t) = m))m)\ 

= 2e- i(B *- Kl) *a*ori*>l^(*)><i|(0i(*)l. 

(13) 

and the reduced density matrix p c (t) of the central sys- 
tem is 

p c (t) - Tr e p(t) 



J2e- i(Ek - El)t a k at (<l>i {t)\4>k (*)>!*) (l\ 



k.i 



(14) 
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The Loschmidt echo L (t) of the central system can be 
calculated as 

L(t) = Tr(p c (t) Po (t)) 

= Tr^e-^^a^ (fr (t) \<f> k (t)) \k) (l\ 
k,l 

xJ^e-^-^K^lm) (n|] 

m,n 

= Tr[^2e-^ E ^a k \ ai \ 2 a: 

k,l,n 

x<<^ (t)\c(> k (t)) t \k) (n|] 

= ^|a fe | 2 h| 2 (^WI^(i)>. (15) 

k,l 

It is clear that if [H c , H ce ] = 0, the decoherence process 
is determined by the initial state of the central system 
{a k } and the time evolution of the {\4>k (t)}}- As shown 
in Eq. (TTTj) . the {\(j) k (t)}} are determined by the initial 
state of the bath (|</>(io))), the eigenstates {\k}} of the 
central system, and the Hamiltonian H ce and H e . The 
eigenvalues {E k } have no effect of the decoherence pro- 
cess. Thus, multiplying H c by a constant does not change 
the L (t) and the diagonal elements of the reduced den- 



sity matrix p c (t) . The time evolution of the absolute 
value of the off-diagonal elements 

\Pc(t) kl \ = \a k at\(<l> l (t)\cl> k (t)), (16) 

is independent of H c . 

Finally, we consider the case that not only [H c , H ce ] = 
but also [H ce ,H e ] = 0. Then, Eq. ([11]) becomes 

\4>k (t)> = e -* M »+ H ^ |0 (t )) = e -Wut e -iH.t | (to)) ) 

(17) 

therefore we have 

(4>i(t)\4> k (t)) = (4>{t Q )\e iH ^e iM ^e- iM ^e- iH ^\cl>{t Q )) 
= ^{t )\e-^ M *- M ^\4>{tv)), (18) 

implying that \p c (t) kl \ and L (t) do not dependent on H e . 
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